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The increasing scale of cryogenic detector arrays for sub-millimeter and millimeter wavelength astro- 
physics has led to the need for large aperture, high index of refraction, low loss, cryogenic refracting 
optics. Silicon with n = 3.4, low loss, and relatively high thermal conductivity is a nearly optimal 
material for these purposes, but requires an antireflection (AR) coating with broad bandwidth, low 
loss, low reflectance, and a matched coefficient of thermal expansion. We present an AR coating for 
curved silicon optics comprised of subwavelength features cut into the lens surface with a custom three 
axis silicon dicing saw. These features constitute a metamaterial that behaves as a simple dielectric 
coating. We have fabricated and coated silicon lenses as large as 33.4 cm in diameter with coatings 
optimized for use between 125-165 GHz. Our design reduces average reflections to a few tenths of 
a percent for angles of incidence up to 30° with low cross-polarization. We describe the design, tol- 
erance, manufacture, and measurements of these coatings and present measurements of the optical 
properties of silicon at millimeter wavelengths at cryogenic and room temperatures. This coating and 
lens fabrication approach is applicable from centimeter to sub-millimeter wavelengths and can be used 
to fabricate coatings with greater than octave bandwidth. 

OCIS codes: 000.2190, 220.0220, 110.0110, 310.1210, 080.2208, 220.3630, 120.4570, 120.4610, 

160.1245, 230.4170 
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1. Introduction 

The development of large format superconducting 
detector arrays for millimeter and sub-millimeter 
astrophysics (e.g. m) has driven the need 
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for high-throughput optical designs that maintain 
diffraction limited performance across the arrays 
(see [2] for a review) . Silicon is an excellent material 
for these applications due to its high index of refrac- 
tion (n = 3.4), low loss-tangent (tan 6 < 7 x 10 -5 ), 
and a relatively high thermal conductivity (k ~ 200 
Wm- l K~ l at 4 I< 0 , which is only a factor of 50 
lesser than high quality OFHC copper [6] and or- 
ders of magnitude higher than plastic). The essen- 
tial development required to realize optical designs 
using silicon optics at millimeter wavelengths is an 
appropriate anti-reflection (AR) coating to mitigate 
the 30% (-5 dB) reflective loss from each optical 
surface. 

AR coatings consist of one or more dielectric lay- 
ers placed on the surface of refractive optical el- 
ements. The index of refraction and thickness of 
the AR layers are chosen such that the reflections 
from the vacuum-AR interface and from the AR- 
substrate interface interfere and cancel. If the op- 
tical path length through a single layer AR coat- 
ing is one quarter wavelength, the refractive index 
of the coating is y/n (where n is the index of the 
substrate), and the coating is free from dielectric 
losses, the cancelation is perfect at one particular 
frequency. Applied to silicon, such quarter-wave 
coatings reduce reflections below 1% over 1.25:1 
bandwidth at normal incidence. 

Wider bandwidths and larger ranges of angles of 
incidence can be accommodated by adding addi- 
tional layers to the AR coating to form a multi- 
layer coating. These coatings require layers with a 
number of different refractive indices. Cryogenic 
applications require that the coefficient of ther- 
mal expansion of the AR layers also be sufficiently 
matched to that of the lens substrate to prevent 
damage upon cooling. The dielectric loss and bire- 
fringence of the AR coating material must also be 
carefully controlled. Locating materials with these 
properties poses a significant challenge to imple- 
menting cryogenic wide band AR coatings for sili- 
con. One avenue to solve this problem is to engineer 
materials with the required dielectric constant by 
cutting sub-wavelength features into the lens sur- 
face mi. Engineered optical materials with prop- 
erties determined by their detailed geometric shape, 
size and orientation are referred to as artificial di- 
electrics or metamaterials. 

Simulated dielectric metamaterial AR coatings 
consist of layers of sub- wavelength holes [|9j, posts, 
or grooves m cut into the substrate to be coated. 
This approach allows the index of the AR coating 
to be precisely tuned by adjusting the geometry of 


the machined features, naturally solves the thermal 
expansion problem, and has loss lower than that 
of the substrate. At cryogenic temperatures sili- 
con has a dielectric loss tangent ~ 100 times lower 
than typical plastics like Cirlex, which have previ- 
ously found use as single layer AR coatings ffH- 
The advantages of silicon metamaterial AR coat- 
ings include precisely controlled indices of refrac- 
tion, inherently matched thermal expansion, and 
significantly reduced dielectric losses. 

This approach has been successfully applied to 
plastic millimeter wave optics through direct ma- 
chining and silicon optics for visible wavelengths 
through patterning of geometric structures which 
perform as a AR coating similar to that found on 
a moth’s eye HM3- At millimeter wavelengths 
Zhang et al. m fabricated an alternative class 
of simulated dielectric AR for silicon composed of 
metal resonant structures supported by a plastic 
substrate that had higher losses than the coating 
described here. Schuster et al. m fabricated sin- 
gle layer silicon micro-machined artificial dielectrics 
and Han et al. m have used silicon immersion 
grating technology to fabricate planar artificial di- 
electrics for THz radiation. Our work is the first 
demonstration of a broad-bandwidth metamaterial 
AR coating on silicon optics with finite curvature 
for millimeter wavelengths. 

The AR coated lenses presented here have been 
developed for the ACTPol [18] project, a polariza- 
tion sensitive receiver for the Atacama Cosmology 
Telescope (ACT, [19]). In Section [2] we present the 
design requirements for the ACTPol optics as they 
apply to this work. The design of the coating is 
described in Section [3] In Section [4] we describe 
the ACTPol lens design. Cryogenic measurements 
of the dielectric properties of silicon samples taken 
from 45 cm diameter ingots at millimeter wave- 
lengths are presented in Section [5] In Section [6] 
we describe the fabrication process of the coating 
on the lenses and present reflection measurements 
of a finished lens. We conclude with a brief dis- 
cussion of the range of applicability of this coating 
technique. 

2. Requirements 

The performance requirements for the AR coated 
silicon lenses described in this work are dictated by 
the scientific goals of the ACTPol project [18] . This 
instrument is designed to detect the faint polar- 
ized signals from the cosmic microwave background 
using three detector arrays fed by independent re- 
imaging optics. Two of the arrays are horn cou- 
pled polarimeters for the 150 GHz band (passband 
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Fig. 1. Left: A ray diagram of one set of ACTPol re-imaging 
optics, which includes three silicon lenses feeding a detector 
array. Right: A histogram of the angles of incidence of rays 
at the surfaces of these three silicon lenses. 

from 125 to 165 GHz), while the third array will 
be a multichroic array including both 150 GHz and 
90 GHz bands (‘90 GHz’ passband from 80 to 110 
GHz) [20] . Extensive design studies demonstrated 
high index of refraction (n > 3 ) lenses were needed 
to achieve high optical quality across the required 
field of view. 

In this paper we focus on the coatings for the 
150 GHz band lenses. These coatings are concep- 
tually similar to those which will be used for the 
broader-band multichroic lenses. The bandwidth 
requirement for the 150 GHz band coatings is to 
have > 99% transmission between 125 and 165 
GHz. The left panel of Figure [I] shows the re- 
imaging optics for a 150 GHz array. The optics 
include three plano-convex silicon lenses with di- 
ameters up to 33.4 cm. The figure shows that rays 
passing through the optical system refract over a 
wide range of angles of incidence. This is quanti- 
fied in the right hand panel of Figure [l] which shows 
that the distribution of angles of incidence is cen- 
tered near 17° and that more than 96% of the rays 
have angles of incidence < 30°. This sets the re- 
quirement that the coatings must be optimized to 
minimize reflections for angles of incidence between 
0° and 30°. 

Since ACTPol is a polarization sensitive exper- 
iment, low cross-polarization is another require- 
ment. Studies of polarization systematics (e.g., 
[ 21 ]) suggest that the CMB temperature to polar- 
ization leakage must be controlled to better than 
1% which corresponds to a requirement that differ- 
ences in the transmission for the two polarizations 
be 0.5%. 

To reduce the thermal emission, all optics are 
cooled to 4K or below. Thus the AR coating must 
be able to withstand cryogenic cycling. 

As described in [j5] high resistivity silicon which 
is available in boules up to 45 cm in diameter, of- 
fers a cryogenic loss tangent tan 5 < 7 x ICR 5 , an 


Fig. 2. The geometry of the two layer metamaterial AR 
coating. The left side shows a cross section through the cen- 
ter of the pillars indicating the design parameters discussed 
in the text. The right side shows an isometric view of the 
structure resulting from making cuts along two perpendicular 
directions. 

index of refraction of n = 3.4, and a relatively high 
thermal conductivity making it an ideal material 
for this optical design. Our approach of directly 
machining metamaterial AR coatings into the lens 
surfaces guarantees that the coatings have low di- 
electric losses and coefficients of thermal expansion 
that are inherently matched to that of the silicon 
lenses. 

3. AR Coating Design 

The AR coating is comprised of two metamaterial 
layers. Each layer consists of an array of square 
pillars on a square grid cut using silicon dicing saw 
blades by making two sets of parallel trenches ro- 
tated 90° relative to each other. It has been shown 
that any electromagnetic structure with an axis of 
n-fold (n>2) rotational symmetry must have po- 
larization independent reflectance and zero cross- 
polarization when light is incident along this axis 
[22] . Hence the 90° rotational symmetry of this 
coating leads to zero cross-polarization at normal 
incidence and low birefringence at oblique incident 
angles. Figure [2] shows a cross-section and isomet- 
ric view of this geometry which is parameterized by 
the depths d \ , cfo , and kerf widths k \ , &2 of the in- 
ner and outer layer respectively and a pitch p. The 
outer layer is cut by using a blade of width &2 and 
the inner layer is cut by making a second set of cuts 
to greater depth with a thinner blade of width k\ . 
We introduce the volume fill factor Vf = (1 — k/p) 2 
to facilitate comparison to published analytic calcu- 
lations based on the second order effective medium 
theory developed by Rytov [23] for the effective di- 
electric constant in the quasi-static limit (p « A) 

[HU2H2S]- 

The nretamaterial layers can be treated as a vol- 
ume distribution of small electromagnetic scatterers 
characterized by electric and magnetic polarizabil- 
ity densities [7]. When the pitch is small compared 
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to shortest wavelength of interest, the fields in the 
layer are homogenous (8j and one can define an ef- 
fective dielectric constant which can be used to pa- 
rameterize the propagation properties of the media. 
This effective dielectric function for the layer, , 
is a function of the density of scatterers in the layer 
[2 7, 28], parametrized by the volume fill factor. As 
the feature size of a composite media approaches a 
significant fraction of the radiation wavelength the 
effective dielectric function becomes frequency de- 
pendent [29] . As the wavelength is further reduced 
the artificial dielectric structure no longer appears 
homogenous and these simple quasi-static consider- 
ations must be augmented to adequately model its 
behavior. In going beyond the quasi-static limit res- 
onant effects, diffraction and scattering can occur 
in the artificial dielectric structure 0 [3DH32] . The 
coatings described here, which operate between the 
quasi-static and diffractive regimes, rely on numeri- 
cal simulations to estimate the relation between the 
volume fill factor and the effective refractive index 
of the metamaterial layer. 

The values of d±, d, 2 , k\, It 2 , and p were cho- 
sen using a three step design process. A prelim- 
inary design was carried out using a classical an- 
alytic model consisting of sheets of dielectric ma- 
terial each with a constant thickness and index of 
refraction [33] . These were varied to minimize the 
reflectance across the band yielding targets for the 
electrical thickness and index of refraction of our 
two layer coating. 

The second step was to translate the index and 
thickness from the simple model into the pitch, kerfs 
and depths for our coating geometry. This step re- 
quired knowledge of the relation between the effec- 
tive index of refraction of an array of square pillars 
and the pitch and kerf of that layer. These relations 
were determined by fitting analytic models to HFSS 
(High Frequency Structure Simulator, [34]) simu- 
lations of single layer coatings with a wide range 
of volume fill factors for a few different choices of 
pitch. The HFSS simulations were carried out us- 
ing Floquet ports and master-slave boundary condi- 
tions to model a box containing a single post of the 
AR coating as an infinite periodic two dimensional 
array of features [34] , 

Figure [3] shows a representative simulation of a 
single layer pillar geometry with insets showing the 
simulated geometry and the simple dielectric layer 
modeled in the fit. The fit constrained the effective 
index n e ff , the effective electrical thickness t e ff and 
the maximum frequency to which our metamaterial 
coating behaves as a simple dielectric. For the pur- 



Fig. 3. A comparison of the reflectance calculated from an 
HFSS numerical simulation of a pillar geometry to fits based 
on a model of a simple dielectric layer. The geometry is 
specified in the inset in the top. In the simulated geometry, 
p = 400 /im, tphys = 220 /im, k = 40 /im and nsi = 3.38. 
A comparison between fits where only the index n is free 
and where both the index n and thickness t are varied is 
shown. The HFFS simulation and the best fit curve differ 
by 5% at 220 GHz (~ 1.36 mm). This corresponds to the 
minimum wavelength, A 5 % (and a corresponding maximum 
frequency fs%) for the specified pitch, p below (above) which 
the coating no longer behaves as a simple dielectric layer. 

poses of our fit we define this maximum frequency 
as the point where the simulated reflectance and 
the analytic fit disagree at more than 5% absolute 
reflectance. We refer to this frequency as / 5 % and 
the corresponding wavelength A 5 %. This indicates 
the transition between specular and diffractive be- 
havior of the coating. 

A second fit which fixed the electrical thickness to 
the physical thickness of the pillars is also shown in 
Figure [3] This fit implies a lower f 5 %. The discrep- 
ancy between the best fit thickness t e ff and the 
physical thickness t p h ys of the metamaterial layer 
arises since the geometric structure of the pillars 
results in fringing of fields at the junction between 
layers causing a perturbative shift between the po- 
sition of the physical interface and the effective lo- 
cation of the junction’s electrical reference plane. 
An analogous effect is encountered in the design 
of metallic waveguide structures 135 , which can be 
analytically treated. 

Figure [4] shows the results from fitting a number 
of HFSS simulations with different pitch and fill fac- 
tor to the analytic model. Plots include the break- 
down wavelength A 5 %, the effective index n e ff, and 
the electrical thickness of a layer with index ^ /nst 
expressed as a ratio of the thickness of a quar- 
ter wave homogenous dielectric layer and a quarter 
wave metamaterial layer. 

The recovered A 5 % is in good agreement with the 
analytic condition derived for a silicon optic in vac- 
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Fig. 4. A summary of results from HFSS simulations of single layers of pillars with varying geometry. Left: The minimum 
wavelength A 5 % for which the single pillar layer is fit by an analytic model to 5% accuracy as a function of pitch p at normal 
incidence. Center: The effective index of refraction n e ff as a function of the volume fill factor Vf obtained by fitting the 
simulated reflectances upto the minimum wavelength X 5 % = 1360 pm for a pitch, p = 400 pm, 40 pm. The line labeled 
“Biber” is a prediction from the analytic Biber model [23] , which corresponds to the quasi-static circuit approximation for the 
media. Right: The thickness t e // for a slab dielectric quarter-wave coating divided by the thickness t p h ys of pillars forming 
an effective quarter wave coating as a function of p/X. The index of the simulated dielectric coating was held fixed at y/nsi 
for these simulations. A quadratic fit is plotted to guide the eye. 


uum m E2 to prevent diffraction from a grating 
array, 

P< 7 , • Q (1) 

[n S i + sin Vi) 

where p is the pitch, A is the wavelength, nsi is the 
index of silicon, and 0j is the angle of incidence. 
This analytic expression shows that the pitch must 
be smaller than implied by the A 5 % calculated at 
normal incidence to minimize diffraction at oblique 
angles of incidence. Consideration of manufactur- 
ing cost and mechanical robustness of the two-layer 
AR coating design favors the largest possible pitch, 
as this choice minimizes the number of cuts required 
to cover a given area and makes the pillars larger 
and therefore stronger. A 435 pm pitch leads to 
acceptable performance up to ~ 175 GHz (~ 1.7 
mm) which is the upper edge of our band. 

Comparing this relation between vj and n e // to 
analytic models }14 . ;24j for the case where the elec- 
tric held of the incident wave is perpendicular to the 
grooves, we found discrepancies that are reduced 
as the pitch decreases and we approach the quasi- 
static limit p/X « 0.1. Given that our coating 
design does not operate in this limit we find the 
analytic models insufficiently accurate for our pur- 
poses. Thus we resort to numerical simulation to 
optimize the geometry of our coating. 

The variation in the electrical thickness compared 
to the physical thickness (Figure [4j right) is shown 
for the case of a quarter wave AR coating with index 
y/nsi ■ Additional simulations show that the elec- 
trical thickness depends on the index of refraction 
of the material on either side of the metamaterial 


dielectric layer. We account for this small effect in 
the final numerical optimization of the multilayer 
coating design. 

With these relations (Figure [4]) in hand we con- 
vert the analytic design of our two layer coating into 
parameters for the pillar geometry. In the final step 
we performed a numerical optimization of the coat- 
ing at 15° angle of incidence using HFSS, with the 
constraint that the cut geometry must match the 
kerf geometry of commercially available dicing saw 
(See Section [6]). This step enables manufactura- 
bility, accounts for any discrepancies between the 
electrical and physical thickness implied by our nu- 
merical simulations, and improves the performance 
at larger angles of incidence. Table [l] gives the 
parameters for the resultant design and Figure [5] 
shows the simulated performance for this coating as 
a function of frequency for several (0°, 15°, and 30°) 


Table 1. The parameters for the ACTPol AR coating. 


parameter 

symbol dimension dimension 



(pm) 

(units of p) 

pitch 

P 

435 

1.0 

kerf* 

ki 

41 

0.094 


k 2 

190 

0.437 

depth 

d\ 

217 

0.499 


d 2 

349 

0.802 

volume fill factor* Vf± 

0.820 



Vf 2 

0.317 



* nominal 

average values 
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simulated performance modeled geometry 



Fig. 5. Left: A simulation of the performance of the AR 
coating designed for the ACTPol lenses at a range of angles of 
incidence. Right: The unit cell geometry (based on measured 
cut profiles) modeled using HFSS. 

angles of incidence. Even at 30° incidence the band 
averaged reflections are at —26 dB and differences 
between the two linear polarization states (TE and 
TM, see [33] for definition) are below 0.5% for all 
frequencies and angles. At 15° incidence, average 
reflections are below —31 dB. The tolerance of the 
design to various possible manufacturing errors are 
described in detail in Appendix A. 

4. ACTPol Lenses 

The ACTPol lenses are cylindrically symmetric 
plano-convex designs in which the convex surface 
is aspheric - a conic section with four perturbing 
terms proportional to the fourth, sixth, eighth, and 
tenth order of the distance from the axis. The 
design optimization procedure for the re-imaging 
optics is similar to that described in [H9| ; the pri- 
mary difference being that the ACTPol optics are 
required to be image-space telecentric to optimize 
the coupling to the flat feedhorn arrays. This was 
accomplished by constraining the chief rays at each 
field point to be near normal incidence at the fo- 
cal plane and by allowing the tilts of lens 2, lens 
3, and the focal plane to vary during the optimiza- 
tion. The resulting design achieves Strehl ratios 
greater than 0.93 across the 150 GHz focal planes 
without accounting for the Gaussian illumination of 
the feedhorns, which effectively improves the image 
quality. All three of the ACTPol optics tubes use 
the same three silicon lens designs with the posi- 
tions and tilts adjusted to optimize the coupling to 
the ACT Gregorian telescope [19| . 

The AR coating approach we have developed con- 
strains the perimeter of the lens designs to enable 
accurate clamping during both machining of the AR 
coatings and cryogenic cycling of the lenses. Fig- 
ure [6] shows a cross-section of the perimeter of an 


ACTPol v28b final 1 50 GHz design - Lens 2 Edge 



Fig. 6. Cross-section of the perimeter of the ACTPol lens 
2, which is designed to enable machining of the AR coating 
and clamping the lens without damaging the AR coating. 

ACTPol lens design. Each lens includes a handling 
ring along its perimeter that is not AR coated and 
is used for mounting the lens during machining as 
well as clamping to cool the lens in the cryogenic 
receiver. The outer corners of the perimeter are 
chamfered to minimize chipping. There are steps 
rising from the perimeter to the lens surface to pro- 
vide clearance between the perimeter clamping re- 
gion and the AR coating. The lens blanks were 
machined by Nu-Tek Precision Optical Corporation 
and achieved 5 /mi tolerances. 

5. Properties of Silicon 

Silicon manufacturing can produce different grades 
of material, such as ultra high purity silicon pro- 
duced by the float zone process [35] that has a neg- 
ligible density of impurities and silicon produced 
by the Czochralski process [35| that has a higher 
level of impurities. As charge carriers and associ- 
ated states introduced by impurities are the cause of 
dielectric loss, use of the highest available purity (as 
inferred from the room temperature resistivity) sil- 
icon minimizes the dielectric losses. Measurements 
of the refractive index, dielectric permittivity and 
loss tangent of various high purity and high resis- 
tivity silicon sample over a range of frequencies and 
temperature have been reported [36f438| . For exam- 
ple, ultra high purity silicon has been measured to 
have a loss tangent of ~ 1 X ICG 5 at room tem- 
perature which would correspond to < 1% loss in 
the ACTPol optical system. Unfortunately, surface 
tension limits the zone-melt purification technique 
used to produce ultra high purity silicon to diame- 
ters below about 200 nun. Therefore, the substrates 
available for large diameter lenses considered in this 
work must be fabricated from Czochralski silicon. 
For our low temperature application we expect the 
bulk conductivity of silicon to freeze out dramati- 
cally reducing the dielectric losses, however, other 
loss mechanisms can persist in the desired design 
band. We validate this general picture of the di- 
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electric loss in silicon by optically characterizing 
the influence of the bulk resistivity as a function 
of frequency and sample temperature. 

Samples produced by the Czochralski process and 
readily available in 450 mm diameter stock with 
bulk resistivities specified to be in the range of 1 to 
> 500 Q-cm were characterized both at room tem- 
perature and 4K using a Bruker 125 high-resolution 
Fourier Transform Spectrometer (FTS) with a Ox- 
ford Cryostat CF continuous liquid helium flow 
sample chamber. The cryostat is equipped with 
pair of 75 pm thick polypropylene windows that en- 
able spectral measurement while allowing the sam- 
ple to be held at a regulated temperature. The sam- 
ples were cut to have a typical thickness of 180 pm. 
double side polished, and placed in a 25 mm di- 
ameter optical test fixture at the focus of an // 6 
beam. The reflections from the two surfaces of the 
sample form a Fabry-Perot resonator for which the 
modeling is relatively simple permitting measure- 
ment of dielectric properties. The silicon samples 
are boron doped (p-type) to adjust the resistivity. 
The > 500 Q-cm resistivity silicon used in the ACT- 
Pol lenses has the minimum dopant level. 

Each sample’s transmission was measured be- 
tween 240 GHz (8 icm) and 18 THz (600 icm) using 
different combinations of sources, beam splitters, 
and detectors for three frequency bands between 
240-450GHz (8-15 icm), 450-2850 GHz (15-95 icm), 
and 2.85-18 THz (95-600 icm). 

The spectral resolution employed, 7.5 GHz (0.25 
icm), fully resolves the sample’s spectral features. 
Sliding stages permit the sample or a reference clear 
aperture to be moved into the FTS beam while in 
the cryostat for in situ calibration. 

The transmission spectrum was modeled as a se- 
ries of homogeneous plane parallel dielectric lay- 
ers (39] . The dielectric function for the silicon was 
approximated by a classical Drude dispersion model 

m- 


where e* = e' r + ie" is a complex function of fre- 
quency oj, the damping rate T, and the contribu- 
tion to the relative permittivity e*^ of higher energy 
transitions. The plasma frequency and the damp- 
ing rates are related by, u 2 p = T / e 0 p, where e 0 is 
the permittivity of free space and p is the sam- 
ple’s bulk resistivity [41]. In this approximation, 
the material’s free carriers are treated as classical 
point charges undergoing random collisions and the 
resulting damping is assumed to be independent of 
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Fig. 7. Measured silicon transmission at room-temperature 
(upper panels) and 4.5 K (lower panels). The left and right 
columns in the panel are for samples specified as 1 < p < 5 B- 
cm and p > 500 fi-cm respectively, the sample thickness is 
~ 180/im thick. The figure contains the measured FTS data 
(red), model (dashed blue), and residual (green). 

the carrier energy. We find this representation suit- 
able to represent the sample’s properties over the 
spectral and temperature ranges of interest. 

The sample thickness is known to ±0.5/rm at 
room temperature and corrected for thermal con- 
traction as a function of temperature [42]. Since 
the fringe rate is proportional to the product of the 
refractive index and sample thickness, uncertainties 
in the sample thickness directly limit the precision 
of the determination of ^(e^). Measurements of 
an optically polished crystalline quartz sample with 
accurately known thickness were used to measure 
a 1% amplitude uncertainty across the entire 
FTS band. This calibration uncertainty leads to 
a corresponding reduction in the measurement’s 
sensitivity to ^(e^). These uncertainties were 
accounted for in fitting the FTS data. These fits 
produced root-mean-square deviations in the range 
of ~0.005-to-0.016 between the model and the 
observation spectra. Representative data for the 
p > 500 and 1 < p < 5 Q-cm samples are shown 
in Figure [7] for sample physical temperatures of 
T p hys = 297 and 4.5 K. 

The results of these fits are shown in Table [2j As 
anticipated the finite plasma frequency and damp- 
ing rate contribute significantly to the room tem- 
perature losses, however, the p > 500 Q-cm sample 
could be approximated by a dielectric constant over 
the spectral range of interest (e.g., T and ui p are con- 
sistent with zero). The influence of free carrier col- 
lisions at room temperature on the optical response 
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Table 2. Silicon Drude Model Fit Parameters. 


p(Tphys 

) Tphys 

°oo 

4 

r/27r 

Wp/27T 

[f2-cm] 

[K] 

[-] 

[-] 

[THz] 

[THz] 

>500 

297 

11.7(1) 

0.0015 

- 

- 


4.5 

11.5(5) 

<0.0008 

- 

- 

l-to-5 

297 

11.6(55) 

0.0046 

0.571 

1.60 


200 

11.5(52) 

0.0049 

0.572 

0.80 


100 

11.4(78) 

0.0028 

0.459 

0.26 


70 

11.4(66) 

0.0011 

0.162 

0.20 


30 

11.4(64) 

<0.0008 

- 

- 


10 

11.4(62) 

<0.0008 

- 

- 


4.5 

11.4(62) 

<0.0008 

- 

- 


becomes more pronounced in the 1 < p < 5 fi-cm 
sample, manifesting as a reduction in transmission 
at low frequencies. However, as the 1 < p < 5 fi-cm 
sample is cooled, bulk conduction is suppressed and 
this effect shifts to lower frequencies than of inter- 
est for millimeter wave applications. The observed 
dielectric parameters are a weak function of tem- 
perature below ~ 30 K as anticipated given expo- 
nential thermal dependance of the bulk resistivity 

m- 

Both samples show ~ 1% shifts in the magnitude 
of the index of refraction (n ~ y/f^o) upon cooling. 
This shift is accounted for in our optical design. At 
low temperature both the high and low resistivity 
sample give an upper limit for the loss tangent of 
tan 5 ~ e^/ej ' x < 7 x 10 -5 . For the ACTPol optics 
at 150 GHz, this corresponds to an absorptive loss 
of < 5% for our optical system comprising three 
lenses. The ACTPol lenses were fabricated from 
p > 500 H-cm silicon as it may reduce loss and 
facilitates room temperature testing. Both high 
and low resistivity silicon would provide acceptable 
performance in 4K cryogenic applications, but 
given the small cost differential high resistivity 
silicon is the natural choice. Although we have 
used high-resistivity silicon, our models and obser- 
vations suggest that the carriers will be frozen out 
at 4K for the range of bulk resistivities considered 
here. We note that the conductivity of silicon 
at room temperature, and therefore the loss, is a 
strong function of UV irradiance. We have not 
investigated this dependence cryogenically. 


6. Coating Fabrication 

The coatings were fabricated using silicon dicing 
saw blades which are available in widths ranging 



licron accurate 
antry . 


dicing saw 


microscope. 


300 nm repeatability 
metrology sensor \ 


Lens Fabrication System 


Fig. 8. The custom, micron accurate, three axis silicon dic- 
ing system used to fabricate coatings on the lenses. The 
labels identify the key components and axes described in the 
Section [6j 


from 10 to hundreds of microns. Mounted on a 
commercial silicon dicing saw machine these blades 
repeatedly cut with micron level precision and cut- 
ting speeds up to several centimeters per second. 
Unfortunately, commercial silicon dicing machines 
are not designed for curved surfaces and cannot ac- 
commodate the large diameter lenses required for 
this project. Therefore we constructed a custom 
three axis silicon dicing system to fabricate these 
coatings. 

The fabrication system is shown in Figure [8j It 
consists of a micron accurate three axis stage on 
which we mount an air bearing dicing spindle and 
a micron accurate depth gauge. The spindle and 
retractable depth gauge are attached to the ver- 
tical stage (Z-axis) which rides on the horizontal 
(Y-axis) stage. A lens is mounted on an aluminum 
mounting plate on a horizontal stage (X-axis) be- 
low the spindle. This plate permits the lens to be 
rotated by 0°, 90°, 180°, and 270° and carries a ref- 
erence wafer that is used in setting up the blades 
prior to cuts. A side looking microscope is mounted 
parallel to the Y-axis to characterize test cuts on 
the reference wafer. Flood cooling water is sprayed 
on the dicing blade while cuts are made to carry 
away debris. A temperature controlled water bath 
is used to regulate the temperature of the spindle 
and maintain the flood coolant and air surrounding 
the dicing system within 1° C. 

Our system does not show any appreciable change 
in cut shape or surface damage when cutting at the 
maximum travel speed of 50 mm/s. We conserva- 
tively operate at 25 mm/s for which it takes a total 
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Lens Front (Convex) 


Reference Cuts (flat side) 


Fig. 9. Left: Photograph of the curved front surface of an 
AR coated ACTPol silicon lens (lens 2) with zoomed view of 
a small patch. Right: Photographs taken with side looking 
metrology microscope of a reference wafer taken prior to and 
after cutting the coating on a lens. 


of 12 hours of machining per lens side. Factoring 
in setup time it takes 6 8-hour days to fabricate a 
single lens. 

Figure [9] shows photographs of one of the ACTPol 
AR coated silicon lenses. The fabricated coating is 
sufficiently robust to permit handling the lenses by 
touching the AR coatings. The manufacture re- 
sulted in less than 10 out of 500,000 posts with 
damage visible by eye. The shape of the cut pro- 
files was evaluated by making cuts on a reference 
silicon wafer before cutting (pre-cut) and after cut- 
ting (post-cut) the lens and measuring them us- 
ing the side looking microscope. Comparison of the 
pre-cut and post-cut measurements (See Figure [9]) 
show that the wide blade cuts repeatably with neg- 
ligible evolution to cut profile while the the narrow 
blade shows some evolution in width in the upper 
third of the cut. Simulations show that this evolu- 
tion, leads to a few tenths of a percent increase in 
the reflectance. The accuracy of the depths of the 
cuts is limited by a 3 micron uncertainty in the zero 
point for the blade depth and a ±7 micron quadru- 
ple warp in the lens. This warp is due to imper- 
fections in the lens mounting plate that stress the 
lens. 

The reflectance of the completed ACTPol lenses 
were measured using a scalar reflectonieter. This 
reflectometer consisted of a tunable narrow band 
continuous wave source which illuminated the flat 
side of lens through a 20° full width half power horn, 
from a distance of 5 cm, and tilted at an angle of 
15° relative to the lens surface normal. A receiver 
consisting of an identical horn coupled to a detector 
diode was placed at the mirror image of the source 
horn relative to the plane defined by the point of 
maximum illumination of the lens and perpendicu- 
lar to a line joining the source and receiver. This 
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Lens Reflectance 15° incidence 


TE polarization measurement 

TM polarization ■ ■ ■ sim. (two sided) 
— sim. (single sided) 


Frequency (GHz) 


-35 


Fig. 10. Comparison between simulations and measurements 
of the fabricated lens. The two sided simulations model the 
AR coating on both sides of a flat silicon sample of thick- 
ness equal to that of the lens at its center. The single sided 
simulations model the AR coating on only one side of a flat 
silicon sample. 


system was calibrated by (1) placing an aluminum 
reflector at the same position as the lens to normal- 
ize the peak reflection to unity, and (2) by remov- 
ing the lens and calibration reflector to measure the 
stray reflections which were found to be negligible. 

The results of this reflectance measurement are 
shown in Figure [TO] for both the TE and TM po- 
larizations. Precise modeling of this measurement 
configuration requires accounting for the interfer- 
ence between reflections from the flat and curved 
lens surface, which is beyond the scope of this work. 
However, the results are expected to be intermedi- 
ate between the reflectance from a single AR coated 
surface (single sided) and the easily simulated re- 
sult for two flat surfaces separated by the central 
thickness of the lens (two sided). These two cases 
are presented in Figure [10] for both linear polar- 
ization states. These simulations incorporate the 
measured cut profiles shown in Figure [9j We em- 
pirically found that moving the lens so that it is 
well centered on the beam increases the interfer- 
ence effects (TM case) while moving the lens off 
center reduces these effects (TE case). Based on 
the reasonable agreement between these simplistic 
simulations and the measurement we are confident 
we are reducing reflections to a few tenths of a per- 
cent over the range of angles of incidence required 
for ACTPol. 

7. Conclusion 

We have described a new approach for AR coating 
silicon lenses over broad bandwidths and a range 
of angles of incidence. Simulations, backed up by 
measurements of an AR coated lens, show that 
the fabricated coating of the lenses presented here 
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can reduce reflections below few tenths of a per- 
cent between 125 and 165 GHz for angles of in- 
cidence between 0 and 30° for cryogenic applica- 
tions. We have developed a micron-accurate 3-axis 
silicon dicing saw facility and are using it to man- 
ufacture AR coated lenses for ACTPol. We have 
also shown that a range of p-type silicon doping 
levels can achieve low loss at cryogenic tempera- 
tures using silicon samples from boules as large as 
45 cm diameter. This approach for implementing 
wide-bandwidth AR coated silicon lenses is appli- 
cable for millimeter and sub-millimeter wavelength 
ranges and can be expanded to wider bandwidth by 
adding additional layers to the AR coating. 
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Appendix A: Tolerances 

The manufacturing tolerances for the eight types 
of machining errors shown in Figure [TT| were evalu- 
ated based on a fiducial design. These include: (1) 
errors in the pitch, Sp; ( 2 ) errors in the centering 
of the deeper groove relative to the shallow groove, 
Sc; (3) errors in the kerf widths, Ski, Ski ; (4) slopes 
in the walls of the grooves parameterized by <$51 
and Ss 2 which represent the distance the upper and 
lower groove tilt inward at the groove bottom if the 
mean width is fixed; (5) errors in the depth of ei- 
ther grooves Sd\, Sdi, ( 6 ) differences in the depth of 
grooves in the two orthogonal directions assuming 
the correct mean depth are parameterized SdLih 
for the upper and SdpR 2 for the lower groove; (7) 
application of this coating to a surface tilted at an 
angle <f> relative to the bottom edge of the saw and 
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Fig. 11. This figure shows the families of machining er- 
rors considered in the tolerance analysis and described in the 
text. The dashed lines show the target shape for the grooves 
and that for the pillars while the white space in the grey 
represents the material as actually cut. Therefore difference 
between the dashed region and the white region represent 
machining errors. For the asymmetric depth and damaged 
row errors the sketches show the geometry of the errors. 
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Fig. 12. The band averaged reflectance at a 15° angle of 
incidence as a function of changes in the pitch. The grey 
region shows range for which the performance is within 3dB 
of the fiducial design. 


where the depth of each groove is measured at the 
center of the groove (e.g., the effect of applying this 
coating at the edge of a curved lens); and (8) rows 
of broken posts. The tolerance to variation in the 
refractive index of the silicon substrate, 5nsi was 
also evaluated. The tolerances for these parameters 
were quantified based on a fiducial model similar to 
what was fabricated but with straight walled pillars 
using HFSS simulations in which the parameters in 
each family were varied separately. For the major- 
ity of these effects the results were distilled to a 
band averaged reflection at 15° angles of incidence. 
The sensitivity to each parameter was quantified as 
the displacement needed to bring about a 3 dB in- 
crease in reflectance. For the case of tilted surface 
and for broken posts the results of a small num- 
ber of simulations were evaluated and compared to 
the fiducial performance. In the remainder of this 
section we discuss these results. 

a. Pitch: Variations in the pitch change the ef- 
fective index of both layers. Figure 12 shows the 
band averaged reflectance as function of variation 
in the pitch. A conservative tolerance of 5p < 3pm 
( dp/p < 0.0075) limits the change in reflectance to 
be less than 3 dB. 

b. Centering: The performance was insensitive 

to errors in the centering of the deeper grooves 
within the shallower grooves. This is consistent 
with both layers of pillars behaving as layers with 
a constant effective index of refraction. 

Kerf: The upper left panel of Figure [13 


c. 


presents the impact of errors in the kerf. Control- 
ling 5k\ and Sk 2 to 3 pm insures a less than 3dB 
degradation in performance. Varying both param- 
eters in the same direction (e.g., 5k\ = dk^) has 
little impact on the performance. We have chosen 
the fiducial design to be near the high end of the 
most favorable region to make the manufacture rel- 
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Fig. 13. This figure shows the impact of errors in the kerf, 
slopes in the sidewalls, overall depth errors, and asymmetric 
differences in depth between the two orthogonal cuts. The 
color scale presents the degradation in reflection at 15° inci- 
dence compared to the simulations in Figure [5] The fiducial 
design is at the coordinate (0,0) in all the plots as highlighted 
by a “+” or arrow. The outer edge of the light blue contour 
(see label) represents 3dB degradation. The horizontal (ver- 
tical) axis represents errors in the inner (outer) layer of the 
coating for the labeled parameters. 


atively immune to blade wear, which could have 
(but did not) narrow the kerf width as machining 
proceeds. 

d. Sloped Walls: The upper right panel of Fig- 

ure [13] presents the impact of slopes in the walls left 
behind by the dicing saw. Conservative estimates 
for the 3 dB tolerance for these two parameters are 
6 pm for <5si and 3 pm for 5s2 ■ Sloping both layers 
by a similar amount results in negligible degrada- 
tion. 

e. Depth: The lower left panel of Figure [13] 

presents the impact of depth errors. The 3 dB tol- 
erances for these parameters are 7.5 pm for 5d± and 
10 pm for Sd 2 . Expressed in terms of optical path 
length the sensitivities are identical. 

/. Asymmetric Depth: The lower right panel of 

Figure [l3]presents the impact of asymmetric depth. 
The 3 dB tolerances for these parameters are 15 pm 
for dd^^ and 12 pm for ddm 2 for asymmetry in 
directions common to both grooves. This error also 
affects the 90° rotational symmetry of the geometry 
resulting in birefringence but the cross-polarization 
can be no larger than the absolute reflection. 
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Fig. 14. Left: The simulated impact of having the surface 
of the lens tilted by 10°. Right: The simulated impact of 
variation in the refractive index of the silicon substrate. 


case produces band averaged reflections below 2% 
which is acceptable though it does produce non- 
trivial cross-polarization. We used simulation to 
estimate that keeping the number of broken posts 
below 1/700 would result in a negligible degrada- 
tion in the overall performance of the coating. 

i. Index: Figure 14, right shows the band 
averaged reflectance as a function of variation in 
the index. A conservative tolerance for the index 
is 5nsi < 0.1 (Snsi/nsi < 0.03). This leads to 
increases in the reflectance of less than 3 dB over 
the fiducial design. This tolerance is weaker than 
that imposed by the optical design. 


g. Application to tilted surfaces: The left panel 

of Figure [141 left compares the performance of this 
coating applied to a surface tilted by 10° compared 
with the performance when applied to a flat surface. 
Fortuitously, this effect improves performance at 
large angles of incidence and has little effect at small 
angles. Applying this coating by making grooves 
parallel to the lens symmetry axis of a curved lens 
will produce acceptable performance. 

h. Groups of broken posts: Assessing the im- 

pact of the broken posts is not straight forward be- 
cause the region and extent of illumination varies 
between the three lenses. The degradation due to 
broken posts would thus be a function of the effec- 
tive aperture area of the lenses and the location of 
the defects in addition to the total area affected. 
Accurately quantifying this is beyond the scope of 
this work. Qualitatively, the degradation would be 
negligible if the fraction of total area affected is 
small and more importantly, the length scale of the 
individual affected areas is smaller than the wave- 
length of incident light. Simulating the impact of 
broken posts is computationally difficult. There- 
fore we resorted to a simplified model consisting of 
a unit-cell containing seven posts in a row with one 
broken off completely. This is equivalent to having 
one seventh of the posts destroyed. This pessimistic 


The 3dB sensitivies derived in this section are 
summarized in Table [3j The band averaged reflec- 
tivity of the fiducial design is -31 dB at 15° inci- 
dence which is significantly better than the require- 
ments for the system. 


Table 3. A summary of the 3 dB sensitivities for all param- 
eters described in the text. 


machining error 

parameter sensitivity 

tolerance 



H m 

(units of p) 

pitch (p) 

Sp 

3 

7.5 x 10" 3 

centering 

6c 

see text 


kerf 

Ski 

3 

7.5 x 10 -3 


Sk 2 

3 

7.5 x 10" 3 

sloped walls 

Ssi 

6 

1.5 x 10" 2 


Ss 2 

3 

7.5 x 10" 3 

depth 

Sdi 

7.5 

1.9 x 10" 2 


Sd 2 

10 

2.5 x 10" 2 

asymmetric depth 

Sd LRl 

15 

3.8 x 10" 2 


ddLR 2 

12 

3 x 10” 2 

tilted surface 


see text 


damaged row 

r 

~ 1/700 

NA 

index 

Sn 

0.1 

NA 



